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1. Introduction, the relevant scales 

With the development of the nanotechnologies, new possibilities recently 
emerged to find electronic properties which reflect the discrete aspect of en- 
ergy spectra. The objects we shall describe in this course are metals which 
will be considered as quantum mechanical objects. Long time ago Kubo 
realized that in a small system of typical size L, the inter- level spacing 
A cx L~'^ may become as large as the temperature. The discreteness of the 
spectrum should manifest itself through a qualitative change in thermody- 
namic properties [ 1, 2]. It turns out that beside A, another characteristic 
energy scale Ec emerges {the Thouless energy), which determines the ther- 
modynamic properties in a crucial way. 

It is the purpose of this course to review the properties of the metallic 
spectra and relate them to physical quantities, transport and thermody- 
namics. In a first approximation, a metal can be considered as a complex 
quantum system which shares universal spectral properties with other so- 
called chaotic systems, like nuclei, molecules, models of billiards. One es- 
sential property is that the energy levels are strongly correlated and present 
the phenomenon of spectral rigidity [ 2]. 

Such a system is described by four length scales: the sample size L, the 
mean free path le which describes the elastic collisions, the Fermi wave 
length \p which depends on the density of electrons and the coherence 
length L^. This last scale is very important because the effects wc aim to 
describe in this course result from the phase coherence of the wave functions 
and thus disappear beyond I/,^ [ 3] . Smaller distances define the mesoscopic 
regime. In this course, we shall mostly consider an electron gas in the 
following limits: 

Xf <^ le <^ L <^ (1.1) 

which correspond to a weakly disordered {Xp <^ le), mesoscopic {L <C L^) 
metal in the diffusive regime {le L). When the disorder becomes so 
large that le is reduced to a length of order Xp, the wave functions become 



5 



6 



G. Montamhaux 



localized on a typical scale ^ called the localization length[ 5, 6]. Here we 
shall study only the diffusive regime where ^ ^ oo. 

In this regime, an electron moves diffusively because it experiences many 
elastic collisions while moving in the sample. The typical distance covered 
by the diffusive particle in a time t varies as r^(t) ~ Dt. The diffusive 
motion is thus characterized by a new time scale td which is the typical time 
for an electron to travel through the sample. It is defined as td = I? jD. 
To this time scale corresponds a new energy scale called the Thouless 
energy: 

where the classical diffusion coefficient D is given by 

D = = -V^ (1.3) 

da 

Te = le/vp is the elastic collision time. ?;f is the Fermi velocity. It is 

interesting to notice that the d.c. residual (T = OK) conductance a can 
be directly related to the Thouless energy. Using the Einstein relation 
a = e^poD, where po is the average density of states, and the Ohm re- 
lation between the conductance G and the conductivity, G = aL'^~^, the 
dimensionless conductance g can be written as 

since the average interlevel distance A is l/{poL'^)- N{Ec) is thus the 
average number of states in a strip of the spectrum of width Ec- We shall 
see that, beside A, the energy scale Ec enters in an essential way to describe 
the spectral properties of mesoscopic systems. Although disordered systems 



quantum ergodic diffusive ballistic 



Fig. 1. Relevant energy scales 



have been mainly studied for their transport properties, there has been in 
the recent years an increasing interest in their spectral properties. This is 
in large part because of the development of nanotechnologies which allows 



Spectral Fluctuations in Disordered Metals 



7 



the fabrication of mesoscopic devices and the detection of small signals. In 
particular, the discovery of persistent currents in mesocopic isolated rings 
has motivated the study of the statistical properties of energy levels in the 
presence of an external parameter, a Aharonov-Bohm (AB) magnetic flux, 
which breaks the Time- Reversal Symmetry (TRS)[ 7, 8, 9, 10]. This 
breakdown of TRS has the effect of changing the spectral properties. It 
has become also possible to measure the spectral correlations by transport 
experiments in which it is possible to do a spectroscopy of the levels [ 11]. 

Another motivation for the study of spectral correlations comes from an 
argument due to Thoulcss which relates the conductance to the typical vari- 
ation of the energy levels in a change of the boundary conditions[ 12, 13]. 
This argument links the transport properties to the spectral properties and 
will be discussed in details. Among the signatures of the phase coherence 
on transport properties, a very important one is the property of Universal 
Conductance Fluctuations. In the mesoscopic regime , the conductance 
is not a self-averaging quantity and its distribution is characterized by a 
universal variance[ 14, 15, 16, 17] 

(G^> - {Gf - C^f (1.5) 

We shall see that this universal variance is closely related to the spectral 
rigidity mentioned above[ 17]. 

The course is organized as follows: the next section presents the mi- 
croscopic model used for disordered metals. It describes non-interacting 
■particles moving in a random potential. It is shown that low energy spec- 
tral properties exhibit the universal features of the Random Matrix Theory 
(RMT) and that deviations exist above the energy scale E^- In section^, a 
semi-classical picture relates the spectral correlations to the description of 
a diffusive particle in a random medium. Section ^ is devoted to the per- 
sistent current in a mesoscopic ring and its relation with the fluctuations 
of the local or global density of states. The curvature distribution and its 
relation with the conductance are analyzed in section ^. Then it is seen 
that the parametric correlations, namely the correlations between energy 
levels at different values of an external parameter like an AB flux, exhibit 
a universal behavior after rescaling of the energy and flux scales (section 
i). 

This course overviews only some aspects of the spectral fluctuations in 
disordered metals, especially those using the semi-classical approach[ 18]. 
The subject has been covered by many aspects on recent review papers [ 
4, 19, 20, 21]. Throughout this course, unless otherwise specified, we 
shall set ?i = l,L'^ = l,A = l and we shall describe spinless electrons. 
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2. Metal as a quantum chaotic system 



The study of spectral statistics in complicated systems has been initiated 
by Wigncr, Dyson, Mehta and others to describe the spectra of nuclei [ 
22, 23, 24, 25, 26]. These authors describe the statistical properties of 
matrix Hamiltonians, with a Gaussian distribution of the elements around a 
zero average. A remarkable feature of the Random Matrix Theory (RMT) 
is that the distribution of levels depends only on the symmetry of the 
Hamiltonian. In particular, if the Hamiltonian is invariant under time- 
reversal symmetry, the statistical ensemble of matrices is invariant under 
orthogonal transformations and is called Gaussian Orthogonal Ensemble 
(GOE). If the system is not invariant under time-reversal symmetry, for 
example in the presence of a magnetic field, the statistical ensemble of 
matrices is invariant under unitary transformations and is called Gaussian 
Unitary Ensemble (GUE). The RMT has been applied to a variety of very 
different physical situations in nuclear, atomic and molecular physics [ 2, 
22, 23, 24, 25, 26]. The remarkable outcome is that it gives a universal 
description of complex spectra and we shall use it to describe the spectra 
of disordered metals. 

The relevance of the RMT for a metal has been first pointed out by 
Gorkov and Eliashberg, in their study of ac response of small metallic 
particles [ 27]. In their work, the complexity was not coming from the 
many-body nature of the Hamiltonian like in nuclear physics but from the 
scattering of electrons on the random shape of the boundaries, the system 
having the structure of a billiard with rough boundaries. Here we describe 
a disordered metal where the electrons are scattered elastically by fixed 
impurities. 

Neglecting the electron-electron interactions, a disordered metal is usu- 
ally described by a one particle Schrodinger equation 

-Ti'^ pA 

nil)= -—{V -i — ftl) + V{Y)i^^Ei^ (2.6) 

2777 Tic 

A is the vector potential. The disorder potential is commonly modeled as: 



{V{v))=Q iyiy)v{y')) = ^±-ll (2.7) 

Te is the elastic collision time. This structure is convenient for analyti- 
cal calculations and the essential physics does not depend on the detailed 
choice of the disorder. The discrete version of this Hamiltonian, called the 
Anderson tight-binding Hamiltonian[ 5], is commonly used for numerical 
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calculations 



tij exp{i— / Adl)ipj 



(2.8) 



where the wave function ipi is defined on the sites i of a lattice. In the sim- 
plest version of this model, the t.y are taken as constant tij — t and V{i) is 
a random variable with a uniform distribution in the range [—W/2, W/2\. 
This Hamiltonian, first described by Anderson, exhibits localization prop- 
erties and a Metal-Insulator transition in 3D [5, 6]. Fig. 2 shows a typical 
spectrum in the diffusive and localized regimes, in the presence of an ex- 
ternal AB fiux (see sect. ||). The goal of this course is to describe the 
statistical properties of such spectra. Let us first recall the usual quanti- 
ties used to measure the level fiuctuations[ 24, 25, 26] (The energies are 
measured in units of the mean level spacing A: 
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Fig. 2. A typical spectrum of a metal in the diffusive regime (left), and in the localized 
regime (right), tp = <j)cj)a is the normalized Aharonov-Bohm flux, see sect. ^ 



— The distribution P{s) of spacing s between consecutive levels. In the 
RMT, it is well described by the Wigner-surmise: 

P(s) cx s^exp(-C;3S^) (2.9) 
where /3 depends on the symmetry of the Hamiltonian. When there is no 
correlation between levels, it has a Poisson behavior: 

P{s) = exp(-s) (2.10) 
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E 



Fig. 3. ^'(•5) and S^(i?) for a metal in the diffusive regime, with (black dots) and without 
(open dots) magnetic flux[ 30]. They are very well described by the Random Matrix 
Theory. 




Fig. 4. Two-point correlation function for a metal in the diffusive regime, with (left) 
and without (right) magnetic flux[ 31]. The continuous line is the RMT prediction A S 
function at the origin is not shown. 



— The two-point correlation function of the Density of States (DOS): 

if(ei,e2) = (p(ei)p(e2))-pg (2.11) 
The average (...) is done on energy. Here the average wih be also made 
on different disorder realizations. In the diffusive regime, these two proce- 
dures are equivalent. If, after averaging, the DOS is constant, the function 
K{€i,e2) depends only on the difference ei — £2- In the RMT and for large 
separation e, it varies as K{e) — l/(/37r^e^). At e = 0, it has a S peak 
which describes the self-correlation of the levels. The full expression is 
simple in the GUE case: K{e) ~ (5(e) — sin^(7re)/(7re)^ (see fig. 4). 

— The number variance Y?(Ei) 

Y?{E) = (SN^iE)) = (N'^iE)) - {N{E)f (2.12) 
measures the fluctuation of the number of levels N{E) in a strip of width 
E. By deflnition, the number variance can be written in terms of this 
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Fig. 5. The structure factor Kit) for a metal in the diffusive regime, with and without 
magnetic fiux. 



two-point correlation function 



^^{E)= / K{ei- €2)deid€2^2 i {E - e)K{e)de (2.13) 
Jo Jo Jo 



In the RMT 



/o Jo 
for E> A 
2 



^ln(i?/A) 



it has a logarithmic behavior: 
Cte 



(2.14) 



Although the number variance is very frequently used in the literature, 
this quantity is not always the most appropriate to describe the correlation 
because it is a double integral of the DOS-DOS correlation function i^r(e). 
Thus the behavior of this quantity at an energy scale E depends on the 
behavior of the two-point correlation function for all energies smaller than 
E. 



The Fourier transform K{t) is called the spectral form factor 
K{t) = -!- / K{e) exp{iet)de 



(2.15) 

The advantage of this quantity is to be directly related to the diffusive 
motion of a classical particle (sect. |^). 

At small times t <^ th where th — 27r?i/A is called the Heisenberg time, 
the form factor varies linearly with time K{t) t/{2ir^l3) and it saturates 
to a constant value l/27r for t ^ oo. This constant is simply the Fourier 
transform of the 5 peak which describes the self-correlation of the levels [ 
28]. Figures 3,4,5 present numerical calculations of the different correlation 
functions, with the Anderson Hamiltonian[ 29, 30, 31]. They are very well 
fitted by the RMT. Using the s upe rsymmetric technique to calculate K{e) 
in the microscopic model (eq. ^.6[ ), Efetov has shown that its expression 
actually coincides with the RMT result [ 32]. 

However, fig. 6 shows that above a given energy which is size and disor- 
der dependent and which has been identified with Ec, the number variance 
exhibits deviations and increases rapidly with energy. These deviations 
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E/A 

Fig. 6. The number variance exhibits deviations to RMT above 

from RMT are not surprising: the Anderson Hamiltonian has a very dif- 
ferent structure from a random matrix. It has many zero matrix elements. 
As we have seen in the introduction, there are other energy scales not in- 
cluded in the RMT. Thus we expect deviations from the universal behavior. 
RMT systems are "ergodic" in the sense that their eigcnfunctions cover uni- 
formly the phase space and have no internal structure. This is clearly not 
the case in a metal at small times or large energies where a diffusive parti- 
cle cannot explore the entire space (fig. 7). More precisely, for small times 
Ie "C td, i.e. large energies E ^ E^, the spatial correlations extend on 
a scale Le = ^DIe = \/h,D/E, so that the system consists of (L/LeY 
independent pieces, fig. 7. Consequently, Y,^{E) ~ {L/LeY {E/EcY/"^ 
[17]. 

In the next chapter, we describe the specific properties of the metallic 




Fig. 7. Schematic difTusion at small times {diffusive regime) and large times (ergodic 
regime) 
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spectra and their deviations from the universal properties given by the 
RMT. 



3. Level correlations in a disordered metal 

We assume that the disorder is weak, so that the energy scales obey the 
hierarchy of eq. |l.l[ Thus the dimensionless conductance g 3> 1. In this 
regime, we want now to account for level correlations within and beyond 
the RMT regime. 

First it is useful to introduce the following important operator known as 
the Green function or resolvent: 

G-^^iE) = (3.16) 

In the site representation, the elements of the Green function are given by 

G^\vy,E) = {r\G^-^{E)\r') = ^ |^l£^lgl (3.17) 

n 

En and 0„ are eigenenergies and eigenstates of the Hamiltonian. The 
Fourier transform G^(r,r',t) = ^ J G^{y,y' ,E)e~^^*dE has the follow- 
ing interesting property: Take a wave function ip{r,t — 0), solution of the 
time dependent Schrodinger equation ( |2.6| ). The time evolution of this 
wave function can be written as: 

V'(r,i)= y"dr'G^(r,r',i)V^(r',0) (3.18) 

So G^{r,r',t) "propagates" a particle from (r',0) to (r, i). It results im- 
mediately that the probability for a particle in r' at t = 0, to be in r at t 
is 

|G«(r,r',i)p = G«(r,r',t)G^(r',r,-i) (3.19) 

If the particle is prepared at the energy ep-, the probability becomes[ 33, 
34]: 

P(r,r',i) = ;^ /p(r,r',w)e-*'^*dw (3.20) 
2tt J 



with 



P(r, r', ^) = ^ {G"ir, r' , + ^)G^{r',r, 6^ - |)) (3.21) 
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The DOS can be expressed in term of the Green function as 

p{E) = -\ j ImG^{r,r,E)dr (3.22) 

It is also useful to define the local DOS as 

p{r, E) = --/mG^(r, r, E) (3.23) 

TT 

The two-point correlation function can thus be written as: 

K{€ue2) = ^j ydrdr'(G«(r,r,ei)G^(r',r',e2))c (3.24) 

where (...)c is the connected part of the average. It is essential to notice 
that the Green function appearing in the DOS and its correlation function 
is a diagonal quantity (in space representation) which depends on the prop- 
agation of particles from some origin to itself. All spectral quantities will 
thus depend only on G(r, r, E). 

Another very important quantity which depends only on the diagonal 
Green function is the return probability to the origin P{t) = P{r,r,t). Its 
Fourier transform P{co) is given by 

P{<^) = ^(G«(r,r,e^ + |)G^(r,r,e^ - |)) (3.25) 

Upon averaging, P{uj) is independent of the position r. 

One sees that K{€i,€2) and P{u)) have similar structures. The goal of 
this section will be to relate these two quantities. To calculate P{t), we 
assume that the metal is well described by a random potential in which the 
electrons experience a difi^usive motion of classical particles so that P{u)) 
is a classical quantity equal to Pci{uj)[ 33, 35]. The probability to diffuse 
from r' to r is thus given by the solution of the classical diffusion equation 
i.e. 

PeKr,r',t) = ^e-^«^*e'q('--'-') (3.26) 

q 

where the diffusion modes q are quantized by the boundary conditions. In 
the limit of an infinite system, it takes the familiar form 

where D is the diffusion coefficient. The return probability and its Fourier 
transform are thus given by[ 36]: 

PcKi) = E^^"^^* ^^'M = E— (3-28) 
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In sect. 3.2, we shall show that the two-point correlation function Kiuj 



£1 — 62) is related to the return probability P{w)- 
3.1. A hriej reminder about weak-localization 

Before describing the structure of the spectral correlations, let us first recall 
a brief qualitative derivation of the first quantum correction to the classical 
conductivity, called weak-localization correction [ 33, 35, 38]. Linear 
response theory shows that the d.c. T = OK average conductivity[ 37] has 
the following structure: 

M = J I drdr'(9,G^(r,r',6f^)9,,G^(r',r,6j.)) (3.29) 

By writing the Green function as a sum of contributions from classical 
paths, as it will be explained in more details in the next subsection (eq. 



3.33), the conductivity has the following structure: 



(a)cx^(B,B*e^(^^-^'')/'') (3.30) 

where Bj and Sj are the amplitude and the action associated to each path 
j. For most of the pairs of trajectories, Sj — Sk > 2nfi, so that their 
contribution to the conductivity cancels in average. The classical conduc- 
tivity is given by the sum of the intensities: Ud oa l^jP- However, 
there is a class of trajectories which can also contribute to the conductiv- 
ity: those which form closed loops. Such a loop can be traveled in clockwise 




Fig. 8. The two closed trajectories which make a closed loop are time-reversal symmet- 
ric. They give an additional contribution to the conductivity called weak-localization 
correction. This correction is proportional to the number of such closed trajectories. 



or anti-clockwise directions. Both trajectories, j and its time-reversed j'^ , 
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have the same action, so that they interfere constructively. As a result, 
in addition to the classical contribution, there is a correction of the form 
Act cx BjB*j, where the sum extends over the closed trajectories. The 

sign of the correction is negative because the trajectories j and have 
opposite momenta. As a result there is an enhanced probability to scatter 
in backwards direction. The conductivity is thus reduced. This effect is 
often called coherent backs cattering or weak-localization. This is a phase 
coherent effect because only trajectories of size smaller than will con- 
tribute to the correction. The amplitude of the correction is proportional 
to the total number of loops of any size. The number of loops of length 
vpt being proportional to the return probability Pci{t), one deduces that 
the total correction is given by: 

^oc-^TPe^Wdi (3.31) 



A^"'^ is the transverse area associated to a semi-classical trajectory. The 
integral is calculated between Xe, the smallest time for diffusion, and r^, 
the time after which the electron looses phase coherence, r^p is given by 
= Dr^p. A magnetic field [ 39] or a Aharonov-Bohm flux[ 40] , by breaking 
the time-reversal symmetry, destroys the weak-localization correction [ 41]. 

In the next subsection, we calculate the spectral function K{uj), by using 
very similar arguments. 

3.2. Semi-classical description of energy levels correlations 

We survey here a method which gives a very physical description of the level 
correlations because it directly relates them to the structure of the classical 
motion of a diffusive particle. More precisely, one can get a relation between 
the spectral form factor K{t) and the classical return probability for a 
particle to return to the origin P{t). This description has been developed 
for diffusive electrons by Argaman, Imry and Smilanskyf 18] from ideas 
originally applied to other various classically chaotic systems[ 42, 43, 44]. 
We present here a brief and slightly different derivation. The first step is 
to write the Green function as a path integral [ 45, 46] 

G«(r,r',t) = ^|"c[x]exp(^S[x]), (3.32) 

where 5[x] = J pdx is the classical action for the path x(t) going from 
x(0) = r' to x(f) = r. The integral extends over all paths going from 
x(0) = r' to x(f) = r. The classical paths given by 5S\x\/5-x. = give the 
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main contribution to the integral (B.32) so that the Green function can be 
written as 

G«(r,r',i?) = ^A,(r,r',i?)e'^^('"-^''^)/\ (3.33) 

3 

This expression results from a stationary phase approximation near classi- 
cal paths[ 47, 48]. The amplitudes Aj result from a Gaussian integration 
around the classical trajectories. Using the definition (3.25) of the return 
probability P{uj), one gets 

P{u;) = ^— (^Aj(r)A^(r)e'[^^(^+t)-5'=(-E-f)]/'i)^ (334) 

This sum is done over closed paths. We have defined Aj{r) = Aj (r,r). 
Moreover, for a closed path, the action is independent of the starting point. 
The main energy dependence is contained in the phase factors and we have 
suppressed it in the amplitudes Aj. It has been argued by Berry that the 
phase factors with large actions such that Sj — Sk ^ ^irh should cancel on 
average[ 28, 49]. As a result, in the absence of time reversal symmetry, only 
the diagonal terms are kept in the sum. Using the relation Tj = dSj/dE 
between the period and the energy dependence of the action for a closed 
path, the phase factors can be expanded as Sj{E + uj/2) = Sj{E) +u!Tj/2, 
so that one obtains: 

Pico) = ^{Y,\A,irre^-^^) (3.35) 
j 

By Fourier transformation, the return probability can be written as: 



P{t)^-^{J2\A,{v)\'Sit-T,)) (3.36) 

3 



After disorder averaging, this quantity is independent of the origin r. We 
now turn to the calculation of the two-point correlation function, that we 
write in a symmetrized form. Using the expression in terms of the Green 
functions, doing the same manipulations as for the return probability, one 
obtains 

k{E,t) = ^ [{p(E+^)p{E-^))e^^*de 



2tt J 2"^^ 2 

3 

3 



j j drdv'A,iv)A*{r')S{t-T,)) 



{} \A,\'Sit~T,)), (3.37) 
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where Aj = J Aj{r)dr. Clearly, the expressions for P{t) and K{t) are 
very similar. To relate them, it is now important to make the following 
distinction. The semi-classical motion is characterized by an ensemble of 
closed trajectories. Pj — \ Aj{r)\'^ is the probability to be on the trajectory 
j with origin r. Many other trajectories follow the same path but they start 
from different points r and r'. This ensemble of trajectories, who have the 
same action, is called an orbit. The probability to be on a given orbit j is 




Fig. 9. a) Two identical trajectories starting from the same origin and belonging to 
the same orbit: these trajectories contribute to P{t). b) Two trajectories with different 
origin belonging to the same orbit: these trajectories contribute to K(t). 



\Aj\'^ = J f drdr' Aj{r)A*{r'). One sees that for one given r, the possible 
choice for r' on the orbit j is proportional to the volume of this orbit. This 
"semi-classical" volume is X^^vpTj were A^~^ is the section of the classical 
orbit and vpTj is its length. Recalling that X^^vp is proportional to the 
inverse DOS, one deduces (this simple argument cannot reproduce here the 
correct prefactor): 

\A,\^^-^P,T, (3.38) 
47rpo 

This result has been obtained for the first time by Hannay and Ozorio de 
Almeida in the ergodic case {Pj = 1) [ 50, 49]. After summation over the 
orbits, comparing eqs. 3.36 and 3.37, one obtains: 

1 

4V2 

This semi-classical expression is believed to be valid for small orbits , 
with period t smaller than the quantum (or Heisenberg) time th = 27r?i/A[ 
18, 49]. Fourier transform of eq. 3.39| gives 

Kic.) = -1^ Im = - J- Rc V , . (3.40) 



k{t) = —tPciit) (3.39) 
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and the number variance is obtained from double integration (eq. 2.12): 



9, . I , Kit) ^,Et, 2 f°° , Pit) ^,Et, 
T?{E) = 8 / dt^ sin2( — ) = - / dt^ sin2(— ) (3.41) 



In the above diagonal approximation, eq. 3.35, we kept the terms where 
k — j. However, in the case where the system is time-reversal invariant, 
for each orbit j there is an orbit which is reversed by time inversion k = j"^ 
and which have the same action. Thus, the non-diagonal terms k = j'^ 
also survive on average. This interference term is not taken into account in 
the classical diffusion so that the semi-classical probability to return to the 
origin is actually twice the classical probability P{t) — 2Pci{t). This de- 
scription is called semi-classicalhecaiise it is essentially a classical picture in 
which the quantum mechanics only appears through a phase coherent clas- 
sical contribution. The two contributions are called in the diagrammatic 
language the diffuson (diagonal) and the cooperon (interference) contribu- 



tions. Eq. 3.40 has been obtained for the first time directly by Altshuler 
and Shklovskii, using a diagrammatic calculation[ 17]. 

In order to take into account the loss of coherence due to dephasing 
events, an exponential damping e~'^* has to be added in the interference 
term of P{t) which cuts the contribution of long orbits, so that the diffusion 



pole in 3.4Cl| has a gap —1^ + Dq^ — > —iu + Dq^ + 7. The inverse scattering 



time 7 is related to the coherence length L^: j = hD/L^ Moreover, the 
semi-classical approximation breaks down for small energy scales uj < A 
so that, even without damping, the divergence in the diffusion pole has to 
be regularized by the transformation uj lo + ijA where 7a is an energy 
scale of the order of A[ 30]. The correlation fimction exhibits clearly two 
distinct regimes: 

When UJ < Ec or T > To, the diffusion is uniform in the sample (fig. 
7): P{t) = 2/P so that K{t) varies linearly in time K{t) = t/(2/37r2). 
By Fourier transform, one has K{uj) = — l//37r^(a; -I- ijA)^ and the number 
variance varies logarithmically: = (2//37r^) ln(£'/7A). This is the ergodic 
regime well described by the R MT. It is obtained by taking only the q = 



contribution in the sum 3.4C . This contribution is called the uniform or 
zero mode. It should be emphasized that this semi-classical approximation 
cannot describe properly time scales close to th or energy scales w ^ A. 

Keeping only the diagonal terms corresponds to the GUE case {(3 = 
2) and the interference terms double the contribution to the correlation 
function (/3 = 1). This explains why the fluctuations in the GOE case are 
about twice as large as in the unitary case. 

In the opposite limit (fig. 7), when 3> i?c or t ^ tjj, the return 
probability depends on the space dimensionality d : P{t) cx V/ {DtY^^ so 
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that k{t) (X t^-d/2ijjd/2 ^^^^ fQj. ^ ^ 



(3.42) 



This leads to a power law dependence of the number variance[ 17] 



It is remarkable that the sign of the correlation function is now dimension 
dependent. In three dimensions, K{lo) varies at large energies positively 
like -\-l/^Juj instead of —l/uj^ for RMT, meaning that the levels tend to 
attract each other at large energies [ 51]. The power law regime of the 
variance has been found numerically recently. However, the deviation of 
the correlation function from RMT is quite hard to observe numerically in 
the metallic rcgimc[ 52, 53, 31] when 5 3> 1. 

Finally, when E > Ti/ts, i.e. t < Tg, the motion is not diffusive anymore 
and becomes ballistic. In this regime, the spectral rigidity is very weakly 
dependent on E[ 54]. 

3.3. Breakdown of time-reversal symmetry, Aharonov-B ohm flux 



The above section has stressed the importance of the interference effects 
in the origin of spectral rigidity. They are affected by the presence of an 
external parameter which breaks the time-reversal symmetry, for example 
a Aharonov-Bohm (AB) magnetic flux: one considers a disordered metal 
having the geometry of a quasi-lD ring of perimeter L pierced by an AB 
flux (f), fig. 10. It is assumed that the vector potential A is a constant (^/L 




(3.43) 




Fig. 10. A ring pierced by a Aharonov-Bohm flux 
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in the ring. The energy spectrum is solution of the Schrodinger equation 
( ^.6[ ) with the periodic boundary condition 'ijj{x + L) = tp{x) where x is the 

coordinate along the ring. A gauge transformation — '4'{x)e^'^ ^ ^'^''^^ 

removes the flux from the Hamiltonian: T-l{(j))'ip = T-l(0)'ip' but the new wave 
function -0' obeys a new boundary condition ^p'{x + L) = ^'{x)e^^'^'^^'^° , 
where (/)o — h/e is the flux quantum. An external Aharonov-Bohm flux is 
thus mathematically equivalent to a change in the boundary condition[ 55]. 
Moreover, it breaks the time reversal symmetry and must lead to a change 
in the spectral statistics, from the orthogonal to the unitary symmetry. 

The transition between the two symmetries for random matrices has been 
considered by Pandey and Mehta [ 56, 57] . They have used an interpolating 
ensemble of matrices of the form H — H{S) + iaH{A) where H{S) and 
H{A) are real symmetric and antisymmetric matrices of dimension N and 
variance u^. {a ~ 0) is the orthogonal case and (a = 1) describes the 
unitary symmetry. These authors have shown that the transition between 
GOE and GUE is driven by the single parameter A = v^o? j t^. For 
Gaussian matrices, the average interlevel spacing is not a constant in the 
spectrum so that the transition depends on the position in the spectrum[ 
56, 57]. For a Gaussian matrix of size N and variance u^, the mean level 
spacing in the band center is given by A = irv /\/N so that the parameter 
which drives the transition is the combination A = No? ji:. For example, 
for small times t ^ th, Pandey and Mehta found that the form factor is 
given by 



and interpolates between GOE (A = 0) and GUE (A ~ 1). To relate these 
parameters to those of our physical problem, let us return to the semi- 
classical description. With the above gauge transformation, the diagonal 
Green function gets an additional phase factor: 



i 

where tp = 0/0o and nij is the winding number of the trajectory j. Fol- 
lowing the same steps as in section one finds 



47rAAt 



(3.44) 




(3.45) 




(3.46) 



and its Fourier transform: 




(3.47) 
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The diagonal terms k = j are unchanged and the terms k — get a 
phase Anrnjif, since, after one turn around the loop, one trajectory picks 
up a phase 2'Kip and its time reversed picks up a phase —2'KLp. Classifying 
the trajectories with respect to their winding number, one finds: 

P(i,</7)=Pe/(i,0)+P„t(t, >f ) 

= -=£= e-™^^V4^*[l + cos(47rm<p)] (3.48) 

^ m— — oo 

The first term is the classical return probability and the second is the 
interference term which oscillates with period 0o/2. We have assumed that 
along the transverse directions, the return probability is time independent 
so that our expression is one-dimensional. The flux dependent part of 
P(t, (^) is shown on figure 11. This figure clearly exhibits the two regimes, 
t < To where the motion is diffusive and the flux dependence is small, and 
t > td where the motion is ergodic and the flux dependence is large. 





Fig. 11. Interference part of the return probability 

Let us first describe the limit where I^td and if -^X so that the sum 
over winding numbers can be replaced by an integral. This is the zero- 
mode approximation which corresponds to the ergodic regime. The return 
probability becomes: 

k{t,^) = -^tP{t,v) = -i^[l + e-i6-^^e^^*] (3.49) 

This result is similar to the one found in the RMT, with the mapping: 
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energy . magEetic flm 



Fig. 12. Variance number T?{E,^p) for different values of the flux i/p and comparison 
with RMT [left). E^(i?, Lp) versus flux for a given energy (right): note that this variation 
is reminiscent of the variation of P{t, ip) vs ip when t ^ tjj, fig. 11 



This mapping is quite natural: Ec/A, which gives the range of the RMT 
correlation, is related to the dimension N of the random matrices and ip 
plays the role of the symmetry-breaking parameter a[ 30]. 

We can then use the exact RMT result to calculate any correlation func- 
tion in a flux for any energy range t < Ec- All the correlation functions are 
thus universal functions of the combination of parameters Ec^'^[ 30, 58]. 
Using a supersymmetric approach, Altland et al. have indeed directly found 
that the two point correlation in the cross-over regime is well described by 
the RMT correlation function [ 59]. 



This dependence ( 3.49 ) of the form factor with the flux has a very sim- 
ple origin: diffusive pairs of orbits accumulate a typical phase (<i>^) = 
((47rTO(/3)^) where {m?) is the typical winding number. In the diffusive 
regime, (m^) = 201/1"^ = 2Eat. (cos$) = e-<*'>/2 = g-ia^^E^v^t ^ The 
symmetry breaking depends on the length of the trajectories. For a tra- 
jectory of time t, TRS is broken for $ ~ 27r, i.e. (pc ~ Xj^JE^t. The 
characteristic flux associated with an energy e is thus 0c ^ \/ ^/ Ec- 

It should be emphasized that this argument holds even for the ballistic 
regime in the non-integrable case, for example for billiards in an AB flux. 
In that case, the motion is ballistic but the accumulation of the phase with 
time is still diffusive[ 60, 61, 62, 63]. 



Using the relation 2.12, we obtain the number variance: 

1 , , 1 . / E^ 



27r^ 7^ 27r^ \ ["f^ + I^Stt'^ Ecip ) J 

This fits perfectly the variation of the number variance with the flux (flg. 
3.3). For -EcV^ ^ A, the number variance is reduced by about a factor 



24 



G. Montambaux 



two, as it is expected from the RMT. 

In order to describe the complete variation of the form factor or of the 
number variance with flux, including the 4>q/2 periodicity, the expression 
for P{t, If) as a discrete sum (3.48) over winding numbers m has to be used[ 
64]. Alternatively, P{t, (p) can also be expanded as a sum over the diffusion 
modes q, (eq. 3.26| ). In the ring geometry, the modes are quantized by 



the periodic boundary conditions[ 64]. For a quasi-lD diffusion, q — 
for the classical term and q = ^{n + 2(p) for the interference term. As a 
result, 



There is an interesting duality between the mode (eq. 3.52 ) and the winding 
number (eq. 3.48) expansions. The zero-mode approximation corresponds 
to a continuous summation over the winding numbers up to infinity, as 
we have done above. On the other hand, the m = winding number 
corresponds to an continuous summation over the modes, i.e. to an infinite 
system where the flux dependence is lost. These two descriptions are related 
by the Poisson summation formula. 

For, E > Ec, i.e. t < td, a diffusive particle has a very small probability 
to encircle the flux. More and more diffusion modes become important 
and the number variance increases with the energy window as a power law. 
However, the flux dependence is still very small. For Ec ^ A, it is found 
to vary as: 

T.^{E,ip) ^ S2(£;,0) - ^ln(l + 4— sin2(27r^)) (3.53) 

TT^ 7A 

One sees that in this regime, the reduction of the number variance is 
much smaller than what could have been naively expected. This is because 
the contribution to the spectral rigidity comes from short time trajectories 
and most of these trajectories cannot enclose a flux. This is clearly seen 
on fig. 11 where the flux dependence of P{t, (p) is small. 

Finally, let us also emphasize that, although the spectrum has the period- 
icity the spectral rigidity has the periodicity </'o/2. The Hamiltonian at 
(j) — (j>o/2 has the GOE symmetry. This is a situation of "false time-reversal 
violation" [ 65] . 

In the next section, we study the physics of the persistent currents which 
is intimately connected to the evolution of the spectral rigidity in a flux. 
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4. Persistent currents 
4..!. Introduction 

Although in the 80's most of the physics of the nicsoscopic systems dealt 
with transport properties, a wide interest recently emerged in the experi- 
mental studies of the equilibrium properties. Among them the search for 
the persistent current of an isolated mesoscopic metallic ring pierced by a 
magnetic Aharonov-Bohm flux (f) = (pcpo . 

Such a ring carries an equilibrium magnetization M which is the deriva- 
tive of the free energy F with respect to the magnetic ficild H: M = 
It corresponds to a current / flowing along the ring, which is given by 



/„ is the Fermi factor associated with the energy c^i. For non- interacting 
electrons, this current is simply the sum of the flux derivatives of each 
individual energy levels. Introducing the flux dependent DOS, p(e, if), the 
current is rewritten as: 



The existence of this current has been predicted long ago by Hund[ 66] 
and first calculated by Bloch and Kulik for a clean strictly ID ring[ 67, 
68]. In 1983, Biittiker, Imry and Landauer proposed that it could be also 
observed in a metallic disordered ring[ 69]. This persistent current has been 
observed for the first time by Levy et al. who measured the magnetization 
of lO'' isolated mesoscopic Cu rings[ 7]. In this experiment, the observed 
magnetization gives only access to the average current of each ring (J) . 

The estimation of this average current raises an interesting problem. 
Prom the spectrum of fig. 2 we see that the energy levels move upwards 
and downwards with the flux in a random way. One can thus expect single 
level current i„ = 5e„/9</) to be zero on average: = and thus (/) = 0. 
In other words, {p{e,(j))) is flux independent. This is because the average 
DOS involves only trajectories of zero length so that no flux is enclosed by 
these trajectories. One then concludes that there is no average persistent 
current (the current is actually exponentially small [ 70, 71]. During the 
last years, there has been a large amount of work to explain the origin of 
a non-zero average current. 

On the other hand, two experiments have been performed on single rings, 
one on Au rings [ 8] , the other on GaAs — GaAlAs rings [ 9]. The relevant 





(4.55) 



26 



G. Montamhaux 



quantity in these cases is the typical current Ityp ~ {{I^) — {I)"^)^^"^ — 

((/2))V2. 

In the foUowing subsection, we present the calculation of the typical and 
average persistent current for diffusive non-interacting electrons. 

4-2. The typical current 

Prom the expression of the current in terms of the DOS, the typical current 
Ityp is given by 

where the correlation function K is the sum of the classical and of the 
interference terms: 

K{u;, if, <p') = Ka{uj, + Ki^ti^, (4.57) 

Similarly: 

P(a;,<^,^') = Pa{uj, ^^) +P,„t(a;, ^±^) (4.58) 

Fourier transforming K{e — e') and using the identity edee*** = —l/t^, 
one obtains straightforwardly: 



i4 

1 r p[ut.v)-p:!im 



87r202 ^3 



dt (4.60) 



where " denotes the second derivative 9^/9i^^. The expansion of the return 
probability in winding numbers gives directly the harmonics decomposition 
of the typical current: 

oo 

^(/2)„sin2(27rmvp) (4.61) 

m=l 

with 
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where P„(t) = -pj^e"'""/*^''*. We use (a ^ 1): 

Jo t» m2"-i 7o 

= 2-i(2a-3)!!^. (4.63) 

where the intermediate integral has the simple meaning of an average di- 
mensionless winding number since w"^ = m'^L'^ /ADt. Taking into account 
the exponential reduction e"'"'* of the return probability at large times 
due to inelastic scattering, one finds that the harmonics are exponentially 
damped by a factor e~"^V'</^<' = ^-mL/L^^ rpj^g harmonics of the typical 
current are then given by: 

, -.O, 96 , En x 9 rw m L , L ,n IT , . „ ,\ 

= ^^fo^ ^' + 2 1; + x^i;) ]^ ' ^ ^'-''^ 

This result should be multiplied by 4 to take the spin degeneracy into ac- 
count, so that the typical current is multiplied by 2. In the limit where 
L L^, the first term in brackets gives the main contribution! 18, 72]. 
It is quite easy to understand qualitatively why the typical current is pro- 
portional to Ec (This proportionality has been first derived in rcf . [ 70] . In 
that calculation, a two-cooperon diagram is missing). It is of the order of 
the charge e divided by the characteristic time which is nothing but td'- 
Ityp oc c/td oc Ec/cfia oc ^jj-j;- It turns out that the calculation of this 
typical current cannot explain the amplitude observed in the experiment 
performed on single Au rings [ 8]. The order of magnitude of the observed 
current is closer to At the moment, there is no theoretical explanation 
for this discrepancy [ 73]. The experiment performed on GaAs — GaAlAs 
rings corresponds to a situation where le — L and the observed current 
corresponds to the theoretical prediction. 

It is also interesting to calculate the single level typical current, ityp. 
Similarly to the typical total current, it is given by 

It differs from the typical total current only by two energy factors. There- 
fore, this implies an additional t"^ term in the integral on time. One finds 
easily: 

(^2)„=8m2(^)2 r^dt (4.66) 
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From the general expression of the integral 4.63, we obtain 

(l2)„, = 87^(^)26-"^/^^ (4.67) 

This expansion in harmonics diverges for — s- oo. This is because the 
semi-classical calculation breaks down for long times t ^ tr, i-e. for large 
winding numbers m ^ a/ A/iJc. It turns out that this expression of the typ- 
ical single level current is only approximate and does not describe properly 
the current resulting from numerical simulations [ 74]. This is because, in 
the calculation of (i^), all energy scales are smaller than A. A calculation 
based on the supersymmetry method is needed [ 75]. 

Finally, it is interesting to notice that a given harmonics of the typical 
total current is g = £^c/A times larger than the single level current. This 
means that although the total current is the sum of the currents of all the 
energy levels in the spectrum, the contribution of the levels nearly cancel 
and only the last E^/ ^ levels contribute to the total current. This is seen 
more precisely with a calculation of the correlation between the harmonics 



of the current taken at different energies: similarly to eq. 4.66, one finds 



^Srn^l^^j Re[e-™V^] (4.68) 

where is the current at energy E. One sees that this correlation function 
decreases on a scale Ec/m'^[ 76]. 

4-3. The average current 

We turn now to the calculation of the average current relevant for the many 
rings experiment [ 7]. 

4-3.1. The canonical current 

The above calculation of the average persistent current, (/) =0, i mpli citly 



assumes that the Fermi level is flux independent, so that in eq. 4.55 , the 
only flux dependence is contained in the DOS. It has been proposed that 
the fact that in each ring, the number of particles N is fixed and not the 
Fermi level, plays an important role and leads to a finite average current [ 
77] . The reason is that the constraint that N — J^^'' p{e, (f>)de is fi xed im plies 



that ep must be flux dependent. This flux dependence in eq. 4.55 being 



correlated to the flux dependence of p(e, (j)) makes the average current non 
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zero [ 78]. This "canonical"0 average can be found by doing an expansion 
around the average Fermi level tp- Denoting by the sample and flux 
dependent chemical potential and ep — {fJ-{(t))) one finds: 

dF, dn, , , 

By definition, the first term of this expression is the grand canonical current 
/p. The first derivative is the number of particles N. As a result: 

dN 

In = le, + ^UAK^) - ^f) (4.71) 

We then use the relation 5^\n ^ — A(57V|^ which expresses that the varia- 
tion of the chemical potential at fixed number of particles is proportional 
to the variation of the number of particles at fixed chemical potential[ 79]. 
After averaging, neglecting (/e^), one deduces the following relation[ 78]: 

{lN)=-^^{5N^{eFA))^-^^fJfjK{e,e\<l^)dede' (4.72) 

Therefore, when the number of electrons in the rings is fixed, the average 
persistent current is finite and is rewritten in terms of the typical sample 
to sample fluctuation in the number of levels below the Fermi energy e^i 
64, 80, 81, 82]. From the above semi-classical analysis, the current can 
be directly written in terms of P{t): 

90 7o 47r2 d(j)jQt 

The integral is known (eq. 4.63) and gives the harmonics expansion 



(In) = 22 sin{Anmip) (4.74) 

m— 1 

with 

(/at)™ = i Ae-^Vi. (4.75) 

TT 00 



* The word canonical is here sUghtly misleading. This calculation is actually performed 
in the grand-canonical ensemble, with a constraint that the average number of particles 
is fixed and flux independent 
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The average current can be reconstructed as[ 83] 

^^n) = — T-TTr — 4.76 

■KcpQ cosnL/Ly — cosiTTip 

In order to take into account the spin of the electrons, the number variance 
should be multiplied by 4 and the interlevel spacing divided by 2 so that 



the current 4.76 is doubled 



This current oscillates with the period 0o/2 and is paramagnetic at small 



flux. This is clearly seen from eq. 4.72: the average current measures the 



change in spectral rigidity when the TRS is broken. The number variance 
decreases when (p is finite so that the current has to be paramagnetic. 
The order of magnitude of this average " canonical" current is smaller than 
observed experimentally for an ensemble of lO"^ Cu rings[ 7]. 

For completeness, we mention here that the magnetism of small disor- 
dered metallic dots - where the external parameter is a magnetic flux in 
the bulk of the dot instead of an AB flux - can be also described along the 
same ideas [ 84]. 



4-3.2. electron- electron interactions 

Another larger contribution to the persistent current may come from inter- 
actions as has been first proposed by Ambegaokar and Eckern, in the frame- 
work of the Hartree-Fock(HF) approximation and using diagrammatic cal- 
culations[ 85, 80] . We present here a simple semi-classical derivation of this 
contribution to the persistent current [ 86]. A different approach, based on 
Density Functional Theory, leads to similar results [ 87]. The HF equations 
read, for each state i 

e,^,(r) =TV'.(r) + F(r)V.(r) + ^ / (7(r - r')|^,(r')| V.(r)dr' 

- E / - r')V*(r')V',(r)V^.(r')dr' (4.77) 

In perturbation to the first order in the interaction parameter [/, the 
shift of the energy levels is thus given by 

e. = e"+E / C/(r-r')|^,(r')nV'.(r)pdrdr' 

-E^-'-. f U{r ~ r')rj{r')^j{r)i^*{r)Mr')drdr' (4.78) 
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where the states V'i are those of the non-interacting system. As a result, 
the total energy Et is now 

Et = E't + \Y.I U{v-v')\^^,{v')\^\^p^{v)\^dvdv' 

^ E*^-.-. / ^("^ - v'Wj{^')^j{^m{^W)dvdv' (4.79) 



where E^ is the total energy in the absence of interaction. The summation 
is on filled energy levels, ai is the spin of a state ipi. In terms of the 
local density n(r) = IVJiCr)!^, the total energy can be rewritten: 

Et = E^ + ]^ j U{r- r')n{r')n{r)drdr' 

E-^-.-^ / - r')V*(r')V',(r)V'*(r)^,(r')drdr' (4.80) 

We now assume a screened Coulomb interaction: U{r — r') = L''^(r — r') 
where U = Aire^/q^p, qtf being the the Thomas-Fermi wave vector. For 
such a local interaction, the Fock term gets the same structure as the 
Hartree term and one has: 

Et = E°+j J n2(r)dr J n^{r)dr (4.81) 

Because of the spin, the Hartree contribution is twice the Fock contribution[ 
85] . The interaction gives a new contribution to the persistent current 

Wc define the local DOS p{r,uj) so that n(r) = 2 p{Y,uj)dijj ( the factor 
2 accounts for spin). The current can be rewritten as: 



9 f 

{le-e) = -Uq^ J wi)p(r, u}2))drdioidw2 (4.83) 

= -^^ J{G''{r,r,u,)G^{r,r,u;2))drdw,du2 (4.84) 

= - — ^ [ P{uJi-u;2)duJidu2 (4.85) 



^_upod_ r 

TT d(j) Jo 



P{t,, 



-dt (4.86) 
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Using again the integral (4.63), one finds easily the harmonics content of 
the HF current: 

oc 

{h-e) = ^(/e-e)msin(47rm(^) (4.87) 

m— 1 

with 

= 16^^[f + m^le-"^/^^ (4.88) 

In the limit L <C L^, the second term in brackets can be neglected and 
one recovers the result obtained for the first time with a diagrammatic 
calculation[ 85, 80, 88]. Starting from eq. 4.86, one can actually recover 
the expression resulting from that calculation 

(/e-e) = ^ . du. (4.89) 

TT d(p ^ Jq 7 + Dq^' — luj 

The quantization of the wave vectors of the diffusion modes in now flux de- 
pende nt q{4 >)- A Poisson summation leads to the above Fourier expansion. 



Eq. 4.88 shows that the persistent current is proportional to the interac- 
tion parameter. It is known that this current is smaller than experimentally 
observed [ 7] . Moreover, inclusion of higher order contributions in the inter- 
action parameter (the so-called Cooper channel renormalization) reduces 
further the amplitude of the estimated current [ 89, 90, 91]. It is one 
order of magnitude smaller than the observed current [ 7]. 



4-.3.3. Local versus global fluctuations of the DOS 

It is interesting to contrast the two expressions for the canonical and 
Hartree-Fock average currents. The canonical current results from a con- 
straint on the conservation of the total number of particles N. From this 
constraint, the current has been rewritten in terms of the flux variation of 
the two-point correlation function of the global density of states. On the 
other hand, the HF average current has been written is terms of the two- 
point correlation function of the localDOS. This must reflect the constraint 
of a local conservation of a number of particles [ 80]. This can be simply 
understood, using a simple argument due to Argaman and Imry[ 87]. 

Suppose that the sample can be divided into pieces i in which the number 
of particles Ni is fixed, due to electrostatic interaction. By definition, 
N = the energies are extensive: F — Fi and fl = ^i- We 

can calculate the current in each box i as we did above for the canonical 
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current. The current can thus be written as 

'.-.=-i:fu.=-i;fu(« (4.90) 

In each box, the local chemical potential /x, has to adjust with the flux to 
keep Ni constant. 

The number of particles in one box is A^, = Then: 

/e-e = /e. +53^|.,(/X^(0)-e^) (4.92) 
i 

Using the relation ^/Ui|jVi = — AJATj]^. , averaging and neglecting {lep), we 
obtain: 

i 

to be contrasted with the canonical current: 

^I^) = -^h,^^^^^'^^^)) ' ^ = E^' (4.94) 

i 

Replacing the discrete sum over the boxes by an integral, one finds (the 
interaction parameter is of order 1): 

= -^± I {n\r)}dr (4.95) 

for the current resulting from a local constraint. This current has exactly 

the structure of the Hartrcc-Fock current (with a prcfactor describing the 
strength of the interaction). On the other hand, the canonical current is 

(^iv) = - 1 (n(r)n(r'))drdr' (4.96) 

and results from a global constraint. 

It is then clear from this viewpoint that the canonical current and the 
Hartree-Fock current have very similar physical origins, one resulting from 
a global conservation of the particle number while the other results from a 
fixed local density due to the electrostatic interaction[ 80]. 
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4-4- Temperature dependence 

Up to now the physical quantities have been estimated at zero tempera- 
ture. The extension to finite T can be easily done by using the identity: 
/ f{^)9{^)de = 2i7rT^^ g{ii^n) where /(e) is the Fermi distribution and 
omegan = (2n + l)7rT. The quantities of interest (except the typical cur- 
rent) have the form 

fdE f dE'T{E-E')f{E)f{E')=An''T^J2J2^{iuJn-ico'J (4.97) 

so that the f-dependent integrands involved in the different expressions of 
the currents should now contain the factor ( ■ T^V^ )^. Since the charac- 

^ smh TzTt ' 

teristic time associated with winding number m is = j Ec, it is clear 
that the temperature dependence of the m*'* harmonics of the current is 
characterized by the energy scale Tm = Ec/ni^. This is consistent with the 
correlation energy which has been found in eq. 4.68, between the harmonics 
of different single level currents. 



5. Conductance and spectrum 

The conductance of a disordered system can be derived from the Kubo 
formula (see below). Its calculation requires the knowledge of both the 
wave functions and the energy levels. We would like to know whether the 
conductance could be also derived only from the knowledge of the spectrum 
and of its correlations. For example, the average DOS does not carry any 
information on the transport or on the degree of localization. But we have 
already seen that the conductance can in principle be obtained from the 
spectral analysis since Ec = gA is the characteristic energy scale at which 
the RMT correlations disappear. We have also seen that Ec = gA drives 
the flux sensitivity of the levels and that a thermodynamic quantity like the 
typical persistent current is proportional to the conductance g. Therefore 
the conductance can somehow be related to the two point-correlation func- 
tion and its variation with an external parameter like an AB flux. A simple 
relation has been first proposed in the 70 's by Thouless, which relates the 
conductance to the curvature of the energy levels. We analyze now this 
relation and more recent developments. 
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5.1. The Thouless formula and its extensions 

A very important step linking the transport to the spectral properties has 
been put forward by Edwards and Thouless [ 12, 13] who argued that the 
electrical conductance can be related to the sensitivity of the energy spec- 
trum to a change in the boundary conditions. Intuitively the more a level 
is localized, the less sensitive it is to the boundary conditions. Therefore 
the conductance should be directly related to a measurement of the sen- 
sitivity to the boundary conditions (we have seen that this is actually the 
case since the typical current is proportional to the conductance) . Thouless 
noticed the very similar structure between the expressions of the conduc- 
tance and of the curvature of energy levels when a change in the boundary 
conditions is introduced ip^x + L) — 'ip{x)e^^ , or, equivalently, an AB flux: 
r] = 2Trip = 2'K(j)/<j)o. On the one hand, the d.c. T = OK conductivity a is 
given by the Kubo formula: 

^ = ^ E M'^i^F - ^o.)5{Ef - ep) (5.99) 

where ea is a single energy level and pafs = {a\p\f3) is the matrix element 
of the momentum operator along the x- direction. On the other hand, the 
curvature of a given energy level a at the origin (77 = 0) is easily found 
from perturbation theory. It is given by: 

In a metallic system, i.e in the presence of moderate disorder, this curvature 
is a random quantity. Thouless assumes first that the matrix elements 
are not correlated with the energy levels e^, , so that the distribution of the 
curvature is roughly the distribution of the l/{ea — £/?). Then assuming 
that the energy levels are not correlated, the distribution of the curvatures 
has the Cauchy form P(c) — (7o/7r)/(7o + c^) with a width 70 given by 

> ^ (5^101) 

On the other hand, assuming again that Pap are decorrelated from the Eq,, 
the Kubo formula gives for the average conductivity: 

^ = ^^(K/5l>o (5.102) 
(...) represents an average over the disorder and the energy levels. Compar- 



ison between the equations 5.101 and 5.102 gives a direct relation between 
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the average dimensionless conductance g — aL'^ /{e^ /K) and the width 
of the distribution of curvatures, known as Thouless relation[ 13]: 



1 To 

2 A 



(5.103) 



The first hypothesis is reasonable and it is correct in the framework 
of RMT. On the other hand, the energy levels are strongly correlated in a 
metal so that the second hypothesis does not hold. Therefore the curvature 
distribution may not have the Cauchy form, see sect. S.l , The conductance 



can also be related to another quantity connected to the flux sensitivity of 
the levels: it has been proposed by Akkermans that the conductance could 
also be related to the averaged square slope of the energy levels instead of 
their curvature[ 92]. The conductance is obtained as 



9d 



A2 



(5.104) 



where i{ri) — — is an average on flux. A similar relation has also 

been proposed by Wilkinson[ 93] and by Simons et al.[ 94, 95, 96] in the 
case where the external parameter, instead of being an AB flux, does not 
break any symmetry, gd measures a global, (^-averaged, property of the 
spectrum, while g given by eq. 5.103 measures a local, tp ^0, property. It 
has been argued that these two quantities which are two different measures 
of the sensitivity of the energy levels to the boundary conditions, should 
actually contain the same information and should be proportional in the 
diffusive regime[ 58]. This aspect will be discussed in section 6.1. 



5.2. Universal Conductance Fluctuations 




a) b) 

Fig. 13. Schematic view of the paths contributing to conductance fluctuations. 



An important signature of the coherent nature of quantum transport is 
the phenomenon of Universal Conductance Fluctuations[ 17, 97]. When 
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a physical parameter is varied, such as the Fermi energy, the magnetic 
field or the disorder configuration, the conductance fiuctuates around its 
average value. These fluctuations are reproducible and are the signature of 
the interference pattern associated to a given impurity configuration. The 
width of the distribution is universal and of the order oi /h[ 97]. As seen 
in section |l|, the origin of this universality is related to the spectral rigidity 
of the spectrum[ 17]: the conductance g is equal to the number of levels 
N{Ec) in a strip of size Ec- The spectral rigidity implies that var[N{Ec)] ^ 
1 so that war [5] ^ 1. This variance can be calculated more precisely directly 
from the Kubo formula[ 17, 97]. The averaged square of the conductance 
contains terms of the form (G^(ri, r;)G^(ri, ri)G^(r2, 4)G^(4, Y2)) (for 
clarity we omitted the gradients). As shown on fig. 13, two contractions 
are possible: y\ = ri,r2 = Y2 and r2 = ri,r2 = y\. The first term a) is 
proportional to J G^(ri, ri, i)G"^(r2, r2, t)(it(iri(ir2. For the same reason as 
in section ri and V2 belong to the same orbit of length vpt. Therefore 
integration on V2 gives a factor proportional to vpt and the corresponding 
contribution to the conductance fluctuation has the form: 

-^cx / tP{t)dt (5.105) 

{gr Jo 

This term has exactly the same structure as the two-point correlation func- 
tion of the DOS. It describes the contribution of the DOS fluctuations of 
the conductance fiuctuations[ 17, 98]. The second term b) is propor- 
tional to / P{r,r' ,t)P{r' ,r,T)dtdTdrdr' . It can be also rewritten in the 
form / tP{t)dt. It describes the contribution of the fiuctuations of the dif- 
fusion coefficient to the conductance fluctuations [ 17, 98]. The integral 
5.105 scales as t|, oc l/(g)^. One then concludes that the fluctuations are 
universal. 



6. Parametric correlations 

6.1. Curvatures distribution 

Thouless assumed that the levels are uncorrelated and found a Cauchy 
distribution for the curvatures. As we know from the RMT, levels are 
actually strongly correlated and repel each other. This repulsion must 
affect the curvature distribution P{c). It is actually easy to find the tail of 
the distribution P{c). When two levels are very close in energy s — > 0, one 
can isolate this pair and treat it in a perturbative way. The distance s(A) 
where A is the perturbation parameter, varies as v? + A^ ^ s + /2s so 
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that for s — > 0, the curvature c varies as 1/s. Since P{s) cx one conchides 
that P{c) — > l/c^+^ for c oo, in contradiction with the Cauchy form[ 
99]. 

The problem of the curvature distribution has been solved recently. In 
the case of pure symmetry, i.e. the perturbation parameter A does not 
break any symmetry ( for example a step potential in the GOE case or an 
AB flux in the unitary case where there is already a magnetic field), the 
curvature distribution is given by: 

Np is a normalization coefficient. This form has first been guessed by 
Zakrzewski and Dclande[ 100] to fit numerical calculations on various model 
exhibiting chaotic spectra: the kicked-top model[ 101], random matrices, 
the kicked rotator and the stadium billiard[ 102], the hydrogen atom in 
a magnetic field. The width 7/3 of the distribution is proportional to the 
average square of the level velocities[ 99, 100]: 



7. = (6.107) 

where i = Afterwards, this curvature distribution has been proven 

analytically by Von Oppen[ 103] for random matrices of the form H(X) = 
H + XK where H and K are random matrices belonging to the same sym- 
metry and A is the perturbation parameter. Recent numerical calculations 
have shown that this distribution is also characteristic of metallic spectra 
when the perturbation parameter is an AB flux (/)[ 74] (see fig. 14). In 
particular, in the limit where (j) —^ 0, the distribution is still the GOE dis- 
tribution (/3 = 1 in eq. 6.106D [ 74]. This has been proven analytically by 



Fyodorov and Sommers who also found that there arc no corrections of 
order l/g[ 104]. 

In the literature, the coefficient 7 entering the Thouless relation is some- 
times taken as the typical curvature. This cannot be the case because, due 
to the 1/c^ tail, the distribution has no second moment. Instead, one can 
choose as a definition of the conductance[ 105, 106, 74, 107]: 

9c ^ ^ (6.108) 



For pure symmetry cases, eq. 3.107 leads to universal relations between gc 
and gd[ 106]: 



GOE del) =71 and 7^ = tt-^^^ ^ gc = Trg^ (6.109) 
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Fig. 14. After rescaling the distribution of curvatures becomes universal. (8x8x8 with 
w = 4,5,6) 



GUE (|c|> = -72 and 72 = 27r-^^^ ^ Ag^ (6.110) 

TT A 

The case of an AB flux lias been studied recently numerically [ 74] and 
analytically with the supersymmetry method [ 104]. The result is that 

gc = 2TTgd (6.111) 

This relation can be simply understood: the distribution of curvatures 
in zero flux is characteristic of the GOE symmetry so that (|c|) is equal to 
the width 7 of the GOE (/3 = 1) distribution (eq. 6.109). But there is a 
finite current only when ip ^ 0, i.e. when the the symmetry has become 
GUE, so that the relation between 7 and {i^ir})) is the relation of the GUE 
(/3 — 2) symmetry (eq. 6.110). Combining these two relations, one gets 
gc = 2TTgd. 

When Lp^ifc, the curvature distribution becomes the one of the unitary 
case with a 1/c* tail, and the second moment now converges. In the cross- 
over regime, it has been found numerically that the curvature distribution 
is quite different than in the pure cases, with a Gaussian tail. The typical 
curvature in the cross-over regime has been found to diverge logarithmically 
at small flux [ 74, 105]: 

(c2(^)) (xg2ln(l/5^2) (g^ll2) 
6.2. Parametric correlations 

We have seen in section ^ that the transition between the GOE and GUE 
symmetries in the diffusive regime is a universal function of the combination 
of parameters i?cV^ or g^"^ . 

Another quantity has been recently introduced by Altshuler et al. to 
characterize the motion of the energy levels [ 108, 94]. It is the auto- 
correlation of the current for a given level. It is defined as C{ip^) = 
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{i{(p)i{(p + f-)), where — is an average on flux. Since the flux drives a 
transition between difl'erent symmetries, the product {i{ip)i{ip + ^p-)) de- 
pends on both if and . It would be translation invariant (independent 
of if) in the case where the parameter Lp does not break any symmetry. 
Here, one has to average on ip so that the correlation function is only a 
function of ip^. One can use the same method as in the previous sections 
to calculate this correlation function. From the definition of the current 
and following the same lines as for the calculation of the typical current, 
one finds [ 108] 

C(^_) = -g^y I K,i{e,e',^)dede' (6.113) 

By Fourier transform and using the relation between K{t) and P{t), one 
deduces[ 109]: 



f2 

1 3^ f P{t,ip_/2) 



2w^ d<pl 



jPit^,, (6.114) 



For < 1, P{t) has the diffusive form: P{t,ip/2) = g-i^^ Ecv^t-it ^ 
that 

^(^-) = g(») (l + ^^^)2 (6-115) 



where b = 7r^C(0) and C(0) = {i^{'p)) = 4£^c/7 is the average single level 
typical current. For large flux ^ <^c = l/\/5' correlation function 
has a universal tail[ 108, 110]: 

C^(^-) = -^ (6-116) 

This form has also been found by Berry and Keating in the case of billiards[ 
111]. It has the same origin, it is related to the diffusive accumulation of 
the phase, even when the motion is ballistic. 

The expression found for C((/?_) does not fit correctly the numerical 
result for small flux. The reason is that the semi-classical calculation holds 
when at least one energy scale is larger than A, which implies here Ec^p^ > 
A. Thus the semi-classical calculation describes only the power-law tail of 
C{ip-), when ip > tpc- On the other hand, it has been show that the small 
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Fig. 15. Universal function C(</3_) 



fliix behavior of C(i^_) has a logarithmic correction: C((^_)/C(0) — 1 cx 
gif \n.l / g<f^[ 105, 74, 75, 112]. This is reminiscent of the logarithmic 
divergence of the typical curvature near the GOE point = 0. 

As the GOE-GUE transition is driven by the unique combination of 
parameters Ecf"^, it is seen that the parametric correlation function 
C{ip-)/C{0) is a universal function of Ecf"^- This universality has been 
stressed by Simons et al. who calculated several other several parametric 
functions, using the supersymmetric method[ 94, 95, 96]. 



7. Conclusion 



In this presentation of the spectral correlations and or their dependence on 
the Aharonov-Bohm flux, we have tried to unify several spectral quantities 
by relating all of them to the same quantity P{t) , the return probability 
for a diffusive particle. Here we summarize schematically the structure of 
several of these quantities, together with the weak-localization correction 
and with the universal conductance fluctuations. This description is based 
on a semi-classical picture in the diagonal approximation and fails when 
all energy scales become smaller than the interlevel spacing. 
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We have recovered in a simple way the main quantities describing the 
persistent currents and the parametric correlations. In these quantities, 
the flux dependence appears through the combination g(p'^, where g is the 
dimensionless conductance. This combination is the typical phase accumu- 
lated by a diffusive particle during the Heisenberg time th- The observed 
order of magnitude of the persistent current in the experiments is still un- 
explained. A promising way could be a better understanding of the role of 
e-e interaction. Several recent works have explored this direction[ 114]. 

In this course, we have restricted ourselves to the description of level 
correlations in the diffusive regime. There has also been a recent activity 
on the study of the correlations at the Metal-Insulator transition which 
occurs in 3D. In the strongly insulating regime, the wave functions of two 
states which are close in energy do not overlap. The statistics becomes 
Poissonian. The Metal-Insulator transition is characterized by a third dis- 
tribution, intermediate between Wigner and Poisson distributions [ 113]. 

Finally, the thermodynamic properties such as the magnetization of bal- 
listic mesoscopic conductors, where the disorder is so weak that 1^ > L 
raises new interesting questions. 
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